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摘要 ”针对 天 环 一 号 (CAT-1: China Astro-Torus 1) 偶 极 场 等 离子 体 物理 实验 要 求 ， 装 置 要 求 中 心 漂浮 
超 导 线 圈 在 无 冷却 和 无 电源 条 件 下 稳定 悬浮 至 少 5hr, 设计 了 超 导 托 举 线圈 与 中 心 漂浮 线圈 耦合 系统 的 悬浮 控 
制 系 统 。 为 确保 1200 kg、5 MA 悬浮 磁体 稳定 悬浮 ， 完 成 了 控制 系统 的 Simulink 模型 建立 和 仿真 ， 基 上 
Routh-Hurwitz 稳定 判 据 ， 研 究 了 PID 控制 策略 对 稳定 控制 影响 ， 确 定 了 稳定 控制 参数 选取 的 范围 : 当 驱 动 托 
举 线圈 的 电源 响应 频率 应 大 于 9 Hz， 输 出 电流 的 误差 小 于 0. 1% 时 ， 理 想 条 件 下 PD 控制 系统 最 有 控制 效果 为 
延迟 时 间 为 0.0463 s、 上 升 时 间 为 0. 1545 s、 峰 值 时 间 为 0. 6283 s、 调 节 时 间 0.0848 s、 超 调 量 5 =1. 6， 
稳 态 误差 err 人 0. 05 mm， 为 电源 相应 频率 大 于 9 Hz， 输 出 误差 小 于 0. 05%。 结 果 表 明 : 基于 PID 能 够 在 较 
短 的 时 间 内 将 悬浮 超 导 环 在 预 设 位 置 ， 且 采用 合适 的 启动 方式 情况 下 ， 可 极 大 的 降低 电路 的 负载 ， 并 为 悬浮 
超 导 悬 浮 偶 的 技术 可 行 性 提供 有 力 支 撑 。 
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Abstract In response to the requirements of the dipole field plasma physics experiment for(CAT-1: China Astro Torus 1), 
the device requires the center floating superconducting coil to be stably levitated for at least $ hours without cooling and power 
supply conditions. A suspension control system for the coupling system of the superconducting lifting coil and the center 
floating coil was designed. To ensure the stable levitation of 1200kg and SMA suspension magnets, the Simulink model of the 
control system was established and simulated. Based on the Routh Hurwitz stability criterion, the influence of PID control 
strategy on stable control was studied, and the range of stable control parameters was determined: when the power response 
frequency of the driving lifting coil should be greater than 9 Hz, and the error of the output current is less than 0.1%, the most 
effective control effect of the PD control system under ideal conditions is a delay time of 0.0463 seconds Rise time is 0.1545s, 
peak time is 0.6283s, adjustment time is 0.0848s, overshoot 6=1.6. The steady-state error err is less than 0.05mm, indicating 
that the corresponding frequency of the power Supply is greater than 9Hz and the output error is less than 0.05%. The results 
show that based on PID, the levitated superconducting ring can be placed in the preset position in a short period of time, and 
under the appropriate starting method, the load of the circuit can be greatly reduced, providing strong support for the technical 
feasibility of the suspended superconducting dipole field device. 


Keywords CAT-1 dipole field device; magnetic levitation; superconducting ring; levitattion control system 
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引言 

从 木星 磁 层 高 效 约 束 空间 等 离子 体 被 发 现 以 来 ，Hasegawa 提出 了 悬浮 偶 极 场 约束 等 离子 体 用 于 聚 等 
离子 体 变 研 究 ， 装 置 具 有 闭合 约束 磁场 。 偶 极 场 等 离子 体 约束 装置 与 其 他 装置 原理 和 结构 的 重大 区 别 是 ， 
如 等 离子 体 约 束 在 悬浮 超 导 磁 体 环 外 侧 ， 呈 稳定 峰值 分 布 ， 悬 浮 体 无 机 械 支撑 ， 无 接触 传 热 ， 高 比 压 (8 ) 
等 特点 。 目 前 国际 上 已 建造 的 悬浮 偶 极 场 等 离子 体 约束 装置 有 2 个 : 一 个 是 MIT 建造 的 LDX”， 另 一 个 是 
东京 大 学 建造 的 RT-1， 两 个 悬浮 偶 极 场 装 置 实验 均 实 验 测 得 在 低频 电磁 扰动 下 等 离子 体 向 内 汕 流 夭 缩 及 
组 织 建立 稳定 性 峰值 分 布 等 ，RT-1 实验 测试 得 局 部 B 比 压 大 于 100%、 电 子 云 约束 时 间 300 秒 所 及 实验 室 测 
到 了 哨 声波 等 。 国 内 悬浮 装置 建造 目前 基本 还 是 空白 ,“ 天 环 一 号 ”(China Astro-Torus No. 1-CAT-1) 悬 
浮 偶 极 场 磁 约束 装置 是 东 华 理 工大 学 团队 国内 首次 提出 设计 研发 ， 用 于 偶 极 场 聚变 高 温 等 离子 体 MHD 稳定 
性 、 输 运 机 制 及 实验 定 标 规律 等 实验 研究 ”。 

超 导 悬 浮 线圈 (Floating Coil-FC) 是 CAT-1 装置 用 于 产生 约束 等 离子 体 的 核心 部 件 ， 超 导 悬 浮 环 重 
量 达 1200 kg 左右 ， 主 要 通过 上 方 托 举 线圈 (Levitation Coil-LC) 电磁 引力 进行 非 接 触 支 撑 ， 超 导 磁 悬 
浮 平衡 与 稳定 控制 对 于 实现 CAT-1 实验 目标 非常 关键 。 根 据 因 绍 定理 ”， 静 电场 、 静 磁场 、 重 力 场 等 不 存 
在 悬浮 稳定 ， 须 引入 非 平方 反比 力 可 以 在 一 定 范 围 内 实现 悬浮 稳定 ， 此 外 工程 上 也 可 以 通过 外 加 反馈 系统 
调节 磁场 参数 以 减弱 不 稳定 性 ”"， 从 而 获得 低 不 稳定 性 悬浮 。 针 对 CAT-1 装置 悬浮 磁体 平衡 与 稳定 性 问 
题 ， 我 们 初步 考虑 的 控制 策略 为 ， 在 水 平和 倾斜 方向 采取 自然 稳定 性 ， 垂 直方 向 采取 外 加 反馈 控制 系统 进 
行动 态 调节 ， 通 过 构建 控制 仿真 模型 来 设计 磁悬浮 控制 系统 的 运行 参数 。 


HH 


1 CAT-1 装置 悬浮 磁体 系统 设计 参数 


CAT-1 装置 物理 目标 是 研究 高 温 高 密 偶 极 场 聚变 等 离子 体 稳定 性 及 消 流 输 运 机 理 ， 以 及 开展 其 等 离子 
体 约 束 性 能 相关 的 物理 实验 研究 ， 装 置 总 体 设 计 参 数目 标 为 等 离子 体 密度 达到 5X 10"m"*， 等 离子 体温 度 达 
到 500 eV， 装 置 真空 室 半径 达 4 m， 超 导 环 外 表面 磁场 过 5T。 总 体 上 领先 国外 现 有 美国 LDX、 日 本 RT-1 
以 及 德国 正在 开发 的 APEX 装置 。 如 图 1 为 CAT-1 装置 3D 结构 及 主机 主要 部 件 示意 图 ，CAT-1 装置 设计 县 
浮 磁体 系统 参数 及 特征 频率 参数 如 表 1 所 示 。 由 托 举 线圈 与 悬浮 磁体 参数 的 设计 可 以 确定 水 平和 倾斜 两 个 
方向 的 处 于 稳定 状态 ， 悬 浮 磁体 应 当 在 预 设 平衡 位 置 的 附近 做 一 定 幅 度 的 振动 ， 允 许 的 振动 区 间 应 当 被 严 
格 限制 。 基 于 CAT-1 悬浮 稳定 总 体 要求 ， 悬 浮 超 导 环 垂 直 振幅 z < 100 mm， 水 平 振幅 er<50 mm， 最 大 倾 
斜 角 a<7/24。 
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表 1 CAT-1 悬浮 磁体 系统 参数 及 特征 频率 参数 


Table.1 CAT-1 Levitation magnet System paramete 


参数 值 
悬浮 磁体 质量 mr 1200.0 kg 
臣 浮 磁体 大 半径 RE 0.5m 
巧 浮 磁 体 小 半径 rr 0.25 m 
巧 浮 高 度 且 2.0 m 
托 举 线圈 半径 RL 1.7 m 
基 浮 磁体 总 电流 IF 5 MA 
托 举 线圈 平衡 电流 IF 3.49 A 
垂直 特征 频率 ff; 0.5 Hz<5 Hz 
水 平 振动 频率 ff。 0.3 Hz< 0.5 
倾斜 振动 频率 f。 0.52 Hz < 1.6 
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图 1 CAT-1 装置 概念 结构 图 


Fig. 1 Schematic diagram of CAT-1 device 


2 PD 悬浮 控制 原理 


1 悬浮 稳定 条 件 
一 般 而 言 ， 无 接触 质点 稳定 存在 某 一 位 置 ， 则 意味 着 质点 沿 任何 方向 产生 位 移 均 有 与 运动 方向 相反 的 
力 使 之 回 到 原 位 置 ， 即 恢复 力 。 对 于 平衡 位 置 悬 浮 线圈 函数 模型 如 为 


CO 089 Ce Xx d7i-me) 


on 


(1) 


式 (1) 中 ,dz 为 托 举 线圈 的 空间 曲线 微 元 ，d 为 悬浮 磁体 的 空间 曲线 微 元 ,了 为 任意 方向 的 单位 矢量 
7 为 悬浮 磁体 上 微 元 指向 托 举 线圈 微 元 的 有 向 矢量 ， 平 衡 位 置 悬浮 线圈 所 受 合 外 力 为 0， 且 mg 为 常数 。 恩 绍 
定理 已 经 指出 上 式 不 能 在 任意 方向 同时 成 立 ， 但 当 苑 为 特定 方向 时 则 可 能 成 立 ， 故 当 悬 浮 线圈 沿 特定 方向 五 
位 移 时 ， 位 移 如 图 2〈a)、'b) 所 示 ， 总 数学 模型 可 以 简化 为 


2 
和， wg $e G3 Xd + $$ x A 0 
人 (2) 
Wo® 人 3 ) dl Xx? x sD 
五 =- ($$ — Xdi+$$ <0 


式 〈2)，z 为 垂直 方 Po er i 以 0F/bz、6RVaer 分 别 表征 悬浮 线圈 在 垂直 和 水 平方 向 
的 稳定 性 。 
除 垂 直 和 水 平方 向 外 ， 对 于 环形 电流 还 需 考虑 倾斜 方向 的 稳定 性 ， 如 图 2(c) 所 示 ， 平 衡 稳定 数学 模 


型 为 
LS 1B (G3) 
式 (3) 中 ， 采 用 归 一 化 的 变量 ， 5 村/94 表 征 县 浮 线 圈 在 倾斜 方 向 的 稳定 性 。 
图 3 所 示 为 超 导 环 基 序 系统 稳定 性 与 位 置 扰动 的 关系 。 图 中 实 线 为 可 稳定 区 域 ， 虚 线 为 不 稳定 区 域 。 
从 图 中 可 见 ， 偶 极 场 悬浮 模型 必须 在 垂直 和 水 平方 向 稳定 性 进行 取舍 ， 而 垂直 方向 为 一 维 ， 水 平方 向 则 为 
二 维稳 定 ， 因 此 悬浮 平衡 点 应 当 在 水 平一 倾斜 稳定 区 域 ， 而 对 垂直 方向 附加 控制 系统 ， 以 保证 偶 极 场 线圈 
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2 悬浮 线圈 运动 方向 分 解 示 意图 


Fig.3 Schematic diagram of the decomposition of the movement direction of the suspension coil 


(a) 垂直 (b) 水 平 (c) 倾斜 


图 3 超 导 环 悬浮 系统 稳定 性 与 位 置 扰动 的 空间 分 量 关系 
Fig.3 the spatial component relationship between the stability and the position disturbance of the superconducting ring 
suspension System 
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2.2 控制 系统 模型 

基于 上 述 CAT-1 的 设计 参数 及 稳定 分 析 ， 本 文 建立 了 垂直 方向 悬浮 磁体 控制 系统 结构 模型 ， 如 图 4 所 
示 ， 其 托 举 线圈 位 于 悬浮 磁体 的 上 方 ， 在 托 举 线圈 中 通 入 适当 电流 对 基 浮 磁体 施加 磁力 以 平衡 重力 。 在 县 
浮 磁体 的 上 下 两 侧 均 有 可 移动 的 限 位 挡 板 以 防 其 失控 损坏 。 本 文 托 举 线圈 磁场 对 悬浮 磁体 的 作用 被 简化 为 
磁体 系统 作为 控制 系统 中 的 一 部 分 来 研究 。 

所 设计 反馈 控制 逻辑 如 图 5 所 示 ， 其 中 参考 电压 由 设计 参数 确定 ,参考 电压 决定 悬浮 磁体 稳定 的 位 置 ; 
控制 系统 采用 经 典 的 PD 控制 ， 传 递 函数 为 QL;， 电源 为 数控 电流 源 以 简化 托 举 线圈 电感 充电 问题 传递 函数 
为 妈 ， 电 源 由 PD 控制 环节 的 数字 输出 控制 ， 将 托 举 线圈 和 悬浮 磁体 的 磁场 作用 整体 简化 为 磁体 系统 以 加 
快 计算 ,传递 函数 为 BB; 悬浮 磁体 的 位 置信 息 通过 激光 位 置 传 感 器 获取 ""， 并 将 位 置信 号 以 电压 形式 反馈 
到 控制 系统 输入 端 ， 传 递 函 数 为 G4。 


PP 
mn 


Ee 国 #3 人 
tts 
DD 支撑 结构 及 基 座 


= 国 Rt 


4 悬浮 控制 系统 结构 


Fig.4 Floating control system structure 
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1) PD 控 
PD 控 


冲模 型 


制 系统 模型 


是 PID 控制 模型 的 简化 模型 ， 相 比 于 PID 控制 ， 响 应 速 


激光 侧 位 
G4 


5 甚 浮 控 制 


Fig. 5 Floating control system block diagram 


数学 模型 如 式 (4) 所 示 ” 


式 (4) 中 , P 为 比例 输 


得 到 PD 控 


2) 电源 系统 模型 


托 举 线圈 处 于 悬浮 磁体 的 强人 磁场 中 ，# 


与 县 


四 节 过 程 中 ， 
型 电流 源 


< 
Ee] 


- oy 


相仿 ' 


3) 磁体 系统 模型 


其中 4 为 电源 的 放大 倍 


磁体 系统 中 ， 


县 


XGA 


P= k,(Vref — Vi) 


d (Vref—Vi) 


D= kg a 


空 制 信号 了 6 


大 
G1= 态 十 Hg= 启 (1 + 5) 


浮 磁体 的 距离 处 于 动态 变化 中 。 为 了 简 
源 ， 影 响 其 输出 的 因素 包括 放大 倍数 、 


者 数 ，w/ 为 电源 响应 频率 。 


平衡 点 附近 可 进行 泰勒 展开 ， 得 到 二 元 拟 合 函 


式 (7) 中 ,zo = 0，I0 = 34943A， 记 AI = I 一 I0、 


在 稳定 状态 下 有 F(ozo) = mg, 经 


mp 


dt? 


电感 不 仅仅 与 上 自 


a 了 | AI+ 守 


AZ=2—2Z0 
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系统 框图 


束 度 更 快 超 调 量 更 小 。 


(4) 


PD 时 域 过 程 


出 , D 为 微分 输出 ，ky 为 比例 放大 系数 ，ka 为 微分 放大 系数 。 将 三 者 县 加 即 可 
制 系统 的 输出 量 gj = P+ D。 在 一 个 稳定 系统 中 ， 


悬浮 磁体 的 位 置 总 是 在 预 设 的 位 置 上 下 波动 最 


后 稳定 在 预 设 位 置 。 再 做 拉 普 拉 斯 变换 ， 获 得 s 域 传递 函数 : 


身 结 构 有 关 ， 
化 控制 过 程 ， 
响应 速度 、 输 出 误差 


数 如 式 〈7) 所 示 
F(Z2) = FUozo) + 0 I) + 


， 可 得 运 


zAz 一 7009 


0(Z — 20) 


云 动 方 程式 (8) 


(5) 


还 与 悬浮 磁体 的 距离 有 关 ， 在 动 


本 文采 用 电流 作为 


空 制 输 出 Li 量 o 


。 电 流 源 的 传递 函数 如 下 ， 


(6) 


浮 人 磁体 仅 在 平衡 点 附近 运动 。 考 虑 到 芒 浮 磁体 的 受 力 与 托 举 线圈 为 线性 


(7) 


(8) 


过 拉 普 拉 斯 变换 ， 得 到 频 域 传 函 式 


关系 ， 在 设计 


核 技 术 
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AZ(s) 
G3= Al7(s) 一 mps? a (9) 
式 (9) 中 ,= 嘲 1。、 z= 叶 s。， 是 在 稳定 性 分 析 中 计算 得 出 。 
4) 激光 测 位 系统 模型 
悬浮 磁体 的 位 置信 息 采 用 激光 测 位 传感器 进行 测量 ， 稳 定 控制 的 可 行 性 很 大 程度 上 取决 于 激光 测 位 传 
感 器 的 准确 度 和 测量 速度 。 激 光 测 位 器 是 基于 光电 学 原理 的 电子 元 器 件 ， 其 电学 特征 主要 输出 幅 值 41/、 响 
应 时 间 及 输出 误差 。 在 本 设计 中 利用 电学 特征 设计 控制 系统 ， 其 传递 函数 模型 为 式 10) 
1 
G4=A 5 
和 (10) 
式 (10) 中 ，4 为 元 件 的 输出 系数 ，w 为 表征 元 件 响应 速度 的 角 频 率 。 
5) 控制 系统 闭环 模型 
基于 上 述 各 控制 环节 的 分 析 建立 闭环 传递 函数 ， 闭 环 传递 函数 为 式 
G1G2G3G4 > (kp + ka* S) (11) 


7 ”1+ 6G16G2G3G4 


式 (11) 中 ， 


a 十 ais 十 azs? 二 a3s3 十 a4S4 


ao = kp 一 Zz0/AiAilo 


al 二 kg 2 /AA olo 


az = mF/AAirlo - zu/ AIAI w1wilo 


as = mr/AA oilo + me/AAio lo 


a4 = mr/AAioiwilo 


式 〈11) 的 零点 为 s 
的 系统 ， 要 求 特 征 方程 仅 在 复 平 国 


上 的 负 半 


平面 有 根 。PD 反 


kp/ka， 特 征 方程 为 ao + ais 二 azs* + a35” + a4S = 0， 对 于 能 
馈 控 制 能 否 实现 悬浮 磁体 的 悬浮 稳定 可 由 


Routh-Hurwitz 稳定 "判别 式 ， 在 理 


an>0 

a a 

a 

aa al 0 

Q4 a2 Q0 >0 
0 Q3 QQ1 


6) 稳定 控制 性 能 指标 


I 系统 的 性 能 
性 能 指标 用 于 描 ; 
次 达到 其 稳 态 
响应 曲线 超过 


生 
Uy 


值 到 达 第 一 个 峰 


间 t。， 即 
值 小 范围 


的 平稳 性 或 系统 的 阻尼 情况 。 稳 态 
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性 能 指标 用 于 描述 抑 种 


症 


基于 稳定 控制 为 前 提 ， 评 价 指标 分 为 动态 怕 
制 系统 从 零 初始 状态 过 渡 到 稳定 状态 的 动态 过 程 。 指 标 包 括 延迟 时 间 tyg， 即 响 
一 半 时 所 需 的 时 间 ; 上 升 时 间 t， 即 啊 应 
值 所 需要 的 时 间 ; 调节 时 
内 所 需 时 间 ， 所 谓 小 范围 通常 是 指 稳 态 值 的 5% 以 内 ;， 超 调 


论 上 存在 稳定 控制 数学 模型 如 式 (12〉 所 示 


(12) 


线 站 
闻 大 ， 即 


FE 能 指标 和 稳 态 和 


次 达到 


仿 值 所 需 时 
| 线 到 达 3 


稳 态 


啊 八 


够 实现 稳定 控 仙 


罗 = ye [14] 
能 指标 两 类 ”。 


问 ， 
保持 在 稳 


一 


峰值 时 


量 


填 之 差 的 百分比 。 如 、 志 、 乌 评价 系统 的 初始 快速 性 ; 妇 描 述 了 系统 


扰动 信号 的 能 力 ， 


0， 即 呵 八 
应 的 总 体 快 ; 
系统 的 收敛 情况 


| 线 的 最 大 偏 


离 量 与 稳 


速 性 ，6 描 述 


了 系统 响 


4。 指标 为 


王 亚 林 ， 刘 云 辉 ， 张 国 书 等 : CAT-1 偶 极 场 超 导 磁 体 的 悬浮 控制 与 仿真 


稳 态 误差 err， 稳 态 误差 的 形成 与 系统 结构 、 参 数 、 外 作用 形成 、 元 件 灵敏 度 、 零 点 漂移 等 各 种 其 他 因素 
有 关 。 控 制 系统 的 性 能 指标 如 图 (6) 所 示 。 


图 6 控制 系统 的 性 能 指标 图 示 


Fig.6 Performance index diagram of control system 
3 Simul ink 控制 仿真 


3.1 CAT-1 磁悬浮 控制 参数 设计 与 仿真 

磁体 系统 环节 的 参数 由 磁体 系统 设计 参数 的 电流 以 及 力学 作用 关系 决定 磁体 系统 物理 参数 如 表 1 所 示 。 
结合 甚 浮 磁体 工作 区 间 的 设计 ， 在 平衡 点 附近 的 土 100 mm 范围 内 求解 受 力 关系 采用 线性 拟 合 对 此 范围 内 的 
受 力 关系 进行 描述 。 在 上 述 讨 论 中 ， 托 举 线圈 的 电流 均 为 总 电流 ， 而 在 控制 过 程 中 ， 电 源 的 输出 电流 为 单 
臣 导 线 内 的 电流 。 托 举 线圈 的 政 数 对 比 于 悬浮 磁体 臣 数 np = 1260 ,假设 为 nj = 1000 进行 计算 ,此 时 
mr1 =327.53 。 为 了 和 迭代 快速 收敛 ， 电 源 参 数 初始 值 参 考 日 本 东京 大 学 RT-1 装置 电源 
(41 =150、wj = 154)。 激光 测 位 环节 的 传递 函数 则 由 产品 的 型 号 有 关 , 参考 日 本 东京 大 学 RT-1 装置 的 激 
光 测 位 传感器 参数 ， 线 性 放大 系数 为 0.05 mm/mV， 响 应 频率 大 于 1 kHz"”， 式 (10) 被 修正 为 G4 = 
20/(1+s/6283.19) V/m。 从 上 述 分 析 ， 可 初步 得 到 CAT-1 控制 模型 参数 ， 如 表 2 所 示 。PD 控制 参数 由 
Routh-Hurwitz 判别 式 求 解 参数 范围 ， 在 此 范围 内 迭代 求 取 最 优 参数 取 值 。 根 据 式 (12) 判别 ， 带 入 表 2 
参数 可 求 得 PD 控制 参数 局 、 和 操 取 值 范 围 为 0.000064 ka< 7.8615 、 
0. 009688< k, <150. 317ka-19. 1205ka-19. 1205ka’。 


表 2 CAT-1 控制 模型 参数 表 
Table 2 Control model 控制 系统 参数 值 parameter table 

1 327.53 

磁体 系统 Zn 9518.92 

Al 20 

激光 测 位 传感器 

油光 测 位 传感器 i 
Ar 150 

| 昌 源 2 


如 图 7，k,。、ka 的 取 值 在 图 7 中 稳定 区 域 ， 取 值 范围 较 大 ， 但 在 动态 控制 过 程 中 甚 浮 磁 体 在 托 举 线圈 
磁场 中 运动 将 会 使 得 磁体 系统 传递 函数 发 生 一 定 的 改变 ， 虽 然 在 CAT-1 装置 垂直 振幅 限制 范围 内 具有 良好 
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的 线性 Sa 但 仍然 会 导致 稳定 区 域 边界 的 移动 ， 因 此 在 实际 的 镶 、 名 应 避免 边界 线 附近 的 值 。PD 控制 输 


出 为 : Vo 


rz Vi)， 其 中 Ver 为 参考 电压 ，V 为 激光 测 位 的 反馈 信号 。 因 此 ， 电 源 受 控 量 为 AL = 工 - I0， 


故 G2 L=150kp(Ve- VD) 。 户 :的 取 值 与 悬浮 磁体 位 置 关系 为 : Ver=IV150kp+20z。 将 
悬浮 磁体 稳定 在 z = 0 处 ， 则 必 有 Vi = 0， 平 衡 电流 为 1 = 34943 A， 电 源 输出 为 Ju = 34.943 A。 


稳定 区 域 


100 
so ] 
0 
0 2 4 6 8 
ka 


图 7 k。、ka 稳 定 取 值 区 域 图 


Fig.7 Stable area with kp、 ka value 


3.2 控制 仿真 模型 


基于 上 述 理论 和 CAT-1 装置 已 经 确定 的 参数 取 值 ， 建 立 Simulink 仿真 控制 模型 "中 。 控 制 过 程 中 ， 参 


考 电 压 , 减 去 激光 测 位 传感器 的 输出 电压 之 后 进入 PD 控制 环节 ， 分 别 作 比 例 放大 和 微分 之 后 求 和 输出 到 


be 


求 和 。 激 》 
指定 了 激光 测 位 环节 初始 输出 ， 系 统 开 关 则 决定 了 控制 系统 的 启动 时 机 。 在 仿真 过 程 中 ， 由 于 电源 必 


须 从 0 输 


Ua 


为 了 更 准 
的 参数 ， 


包 源 的 信号 输入 端 ， 用 于 调节 电源 的 输出 。 电 流 输出 与 初始 电流 仅 加 之 后 进入 磁体 环节 ， 磁 体 环 节 的 传递 
函数 是 基于 初始 平衡 状态 假设 基础 上 求 得 ， 因 此 其 
EC 测 位 环节 的 输出 作为 控制 系统 的 反馈 回 到 PD 控制 环节 的 入 口 ， 完 成 团 环 。 


BE 流 输入 需要 与 初始 电流 求 差 ， 其 输出 需要 与 初始 位 置 


出 开始 上 升 电流 ， 这 将 引起 仿真 的 初始 状态 不 稳定 。 在 真实 的 CAT-1 悬浮 控制 过 程 中 ， 
悬浮 磁体 由 机 械 结构 送 至 指定 位 置 ， 托 举 线圈 电流 达到 平衡 电流 后 ， 机 械 结 构 释 放 悬 浮 磁 体 。 因 此 ， 


确 地 对 控 制 系统 进行 仿真 ， 系 统 开 关 的 设置 是 必要 的 。 如 图 8 中 具有 红色 边缘 的 模块 包含 未 确定 


这 些 参数 的 值 将 在 仿真 过 程 中 由 matlab 程序 传 入 。 


327.53 
[ni 
95 
To lg a 12007 ~9518.92 


图 8 Simulink 控制 仿真 模型 


Fig. 8 Simulink control model 


ka 参数 与 控制 系统 的 性 能 
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在 PD 控制 过 程 中 ，P 分 量 与 悬浮 磁体 在 托 举 线圈 
直接 表征 悬浮 磁体 的 位 置 ; D 分 量 为 巧 浮 磁 
立 置 调 节 过 快 而 导致 速度 失控 或 位 置 反 冲 ， 了 
难以 收敛 ， 取 值 过 小 则 可 能 发 生 速度 失控 或 反 ? 
悬浮 磁体 的 快速 收 公 。 为 得 到 合适 的 局 、Fa 参 数值 ， 簿 
区 域内 进行 扫描 ， 获 得 kK,、ka 与 控制 系统 的 性 能 关系 。 虽然 k,、ka 参 数 稳定 取 值 
关系 统 中 难以 支持 控制 
的 5V 工作 限 1 


里 


参数 部 分 稳 


但 在 实际 的 探 外 
遍 离 幅度 和 弱 : 


浮 磁 体 人 


大 会 导致 晤 浮 磁 
单 的 ，Ku 取 值 过 大 将 导致 悬浮 磁体 移动 速度 过 慢 
此 ， 适 当 的 操 、fa 取 值 有 利 卫 


人 


9 


体 


超 高 电压 信号 的 输出 ， 


磁场 中 的 位 置 (z) 线性 相关 P = ky (Vref 一 412)， 


本位 置 在 时 间 上 的 微分 即 惹 浮 磁体 的 速度 。 若 k, 取 值 
区 值 过 小 则 会 使 得 控制 系统 难 


时 高 ! 


同 


庆 ， 限 定名 的 最 大 取 值 范围 


BB 压 输 入 也 会 损坏 1 
为 : kpmax=5/ (Zmax—Zmin) Al, 其 中 ziiox 为 


过 
同 


以 收敛 。 
:过 大 。 
8 模型 ， 对 其 
区 域 非常 大 ， 
此 ， 根 据 悬 


娠 . 


XGA 


多 


I 用 


电源 。 因 


浮 磁体 位 置 的 最 大 取 值 ， 根 据 经 验 被 限定 为 0. 01 m; zi 为 悬浮 磁体 位 置 的 最 小 取 值 ， 被 限定 为 -0. 01 m; 
4 为 激光 测 位 传感器 的 放大 倍数 , 取 值 为 20v/m, 故 kymax=1. 25。 为 得 到 合适 的 k,、ka 参 数值 , 在 0.000064 < 


k, < 1.25、0.00097 < ka < 7.8615 内 ， 利 用 图 


指标 ， 性 能 指标 如 图 9 所 示 。 


时 间 t (s) 


| 
一 一 一 一 


?sb 
147 5 
4.28571 六 一 bp 
嘴 
12.85714 怕 t 
最 史 
入 
1142857 < ”一 人 
吧 
加 -一 0 
0.05 01 0.15 0.2 山 
ky pl 
图 9 ka 值 与 控制 系统 各 性 能 指标 关系 


8 模型 ， 对 参数 进行 扫描 ， 获 得 k,、ka 与 控制 系统 的 性 能 


Fig.9 Relationship between kavalue and performance indexes of control system 


扫 


过 程 
将 稳定 位 置 调节 为 z = 0.005 m; 时 间 上 = 10s 时 , 仿真 结束 。 为 了 加 快 仿真 过 程 ， 
时 间 段 判断 系统 是 否 收敛 。 对 于 在 t = 3~8 s 时 间 内 未 能 收敛 的 仿真 过 程 将 被 忽 


1, 在 仿真 时 间 t = 2 s 时 打开 系统 开关 ;时 间 t= 3s 时 以 阶 跃 信号 的 形式 变换 参考 1 


过 分 长 收敛 的 时 间 没 有 应 用 的 意义 。 参 数 扫描 初步 结果 为 : 


(1) 在 0.000064 < kk, < 1.25、0.00097 < ka < 7.8615 内 的 大 部 分 
章 。 当 ka 大 于 2.43707 时 ， 无 论 所 在 容许 范围 


图 中 左上 角 


区 域 可 以 实现 稳定 控 人 


系统 在 5 s 达到 稳定 收敛 。 


(2) 对 
的 值 变 大 时 ,位置 控制 逐渐 演变 为 速度 控 外 
使 其 难以 达到 稳定 状态 ， 而 局 在 扫描 的 区 域内 
电源 输入 端 有 尽 可 能 好 的 分 辨 能 力 ， 馈 的 值 应 尽 可 能 的 大 ， 因 此 令 局 = 1.25，ka 的 值 通 
电压 变化 将 引起 
j 当 ka 太 大 时 控制 过 程 将 会 剧烈 拌 动 ， 难 以 在 指定 时 间 内 实现 稳定 甚 浮 。 
取 需 要 从 动态 性 能 指标 和 稳 态 性 能 指标 来 考虑 ，kg 值 与 各 指标 关系 如 
随 着 ka 值 的 增 大 ，tg、 刀 、t、ts 均 有 不 同 程 度 的 增 大 ，ka 大 于 0.073 之 后 ， 无 超 调 


(3) 为 了 使 ! 
过 参数 扫描 的 方式 获得 ， 当 kg 过 小 时 ， 参 考 晶 
的 挑战 ， 
ka 的 最 优 值 选 ] 
系 表明 ， 
0.4 之 后 ， 控 制 过 程 不 用 


| 0.000064 < kK, < 1.25、0.00097 < ka < 2.43707 区 域 进 行 更 为 精 名 


三 是 
总 征用 


电压 广 e 刀 
其 中 在 时 间 t= 8~10s 


区 域 均 不 


略 ， 因 为 在 实际 的 工程 中 


ob nr 


能 实现 稳定 控制 ， 仅 有 
内 如 何 取 值 ， 均 不 能 使 控制 


二 央 
里 ， 


巨大 的 超 调 


和 有 峰值 ， ka 大 于 2.0 之 后 超 调 量 与 稳 态 误差 重合 ， 
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的 扫描 ， 结 果 表 明 ， 当 ka 


关 ， 控 制 系统 对 巧 浮 磁 体 的 速度 非常 敏感 ， 而 造成 悬浮 磁体 抖动 ， 
够 使 系统 收敛 。 


相应 地 对 电源 的 输出 也 是 巨大 


名 


9 所 示 。 图 中 关 
量 ， 当 ks 大 于 
即 其 泽 磁 体 在 稳 态 值 附近 振动 ; 


E。 
= 二 


核 技 术 


20XX, XX: XXXXXX 


Ku 大 于 2.2 之 后 ， 控 种 


3.4 磁悬浮 控制 仿真 


由 


站 


基于 上 述 控制 参数 的 结论 ， 在 1000 s 的 控 


悬浮 实验 中 的 超 调 量 处 于 同上 


4 总 结 


本 文 将 托 举 线圈 


仿真 过 程 初步 结果 如 表 3 所 示 。 


和 巧 浮 磁体 作为 整体 建立 其 在 控制 系统 


省。 


判 仿真 中 ， 忽 略 了 | 


| 过程 不 能 在 仿真 时 间 内 达到 稳定 。 综 上 所 述 ，ka 最 优 值 为 0.073。 


电源 输出 的 纹 波 有 影响， 控制 超 调 量 与 RT-1 
里 想 情 况 基本 与 RT-1 的 悬浮 控制 具有 相同 的 效果 。CATL1 装置 的 控制 


的 传 函 模型 , 建立 了 闭环 PD 控制 系统 ,分析 


PD 控制 参数 与 系统 稳定 的 关系 。 通 过 参数 扫描 获得 最 优 PD 控制 参数 。 得 到 结论 如 下 : PD 控制 系统 在 合适 
的 电源 响应 频率 应 大 于 9 Hz ， 输 出 电流 


的 参数 设置 下 能 够 实现 CAT-1 磁 基 浮 的 稳定 控制 。 驱 动 托 举 线圈 


的 误差 应 满足 式 4 6 才能 保证 悬浮 磁体 的 稳定 悬浮 ， 进 一 步 提高 控制 稳定 


0. 05% 水 平 才能 与 RT-1 控 


1.69%。 
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